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INTRODUCTION 


One of the challenges in the development of a hydrogen-fueled, combustion ramjet 
engine is to assure that combustion is maintained over the entire pressure- 
temperature regime encountered* A major problem has been the uncertainty in ignition 
delay times for the candidate fuel. To cope with this problem, the use of combustion 
initiators is being seriously considered. The addition of a combustion initiator 
would ensure combustion during critical regions of pressure and temperature, which 
are dependent on vehicular speed, altitude, and other flight characteristics. One of 
the most promising combustion initiators being considered presently is silane, SiH^. 

Experiments in the Langley Research Center Hypersonic Propulsion Test Cells and 
at contractor facilities have demonstrated that the addition of silane to hydrogen 
fuel enhances combustion and also helps sustain stable combustion (ref. 1). A chem- 
ical kinetic mechanism for the silane-hydrogen system was suggested by Beach et al. 
(ref. 2) • Calculations with this mechanism also showed that significant reduction in 
ignition delay times is possible when silane is added to hydrogen-air mixtures. 

The objective of the study reported here was to quantitatively characterize the 
ignition characteristics of silane-hydrogen mixtures and to use the results to test 
the mechanism proposed by Beach et al. (ref. 2). Tie experimental study was carried 
out in a chemical kinetic shock tube. Ignition delay times were measured behind 
reflected shock waves for two silane-hydrogen-oxygen-nitrogen mixtures at test con- 
ditions which duplicate pressures and temperatures which are likely to occur within a 
supersonic ramjet combustor. Tie measured ignition delay times were compared with 
the ignition delay times predicted by the si lane -hydrogen reaction model. The the- 
oretical comparisons were made with a computer program which simulates the chemical 
kinetics behind reflected shock waves (ref. 3) • The thermochemical data used in mak- 
ing the theoretical comparisons were taken from the literature, primarily JANAF 
(ref. 4) • For some fragmentary compounds deemed necessary to describe the silane 
combustion process, thermochemical properties were estimated. 


EXPERIMENTAL APPARATUS AND MEASUREMENTS 

All the experimental measurements were made behind reflected shock waves in a 
stainless steel tube with an inside diameter of 8.9 cm. The test section was 671 cm 
long. The present arrangement differed from that shown in references 5 and 6 in that 
the dump tank was removed and replaced with a flat stainless steel plate. Fast- 
response, transient pressure transducers were spaced every 30.48 cm along the last 
122 cm of the test section. Hie test section and instrumentation locations are shown 
in figure 1 • The flat plate located at the end wall of the shock tube was machined 
to accept a transient pressure transducer at the axial centerline of the shock tube. 
The pressure transducer placed here actually protruded 3 mm into the approaching 
shock wave and incoming flow. Ignition measurements were made at this location. 
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Figure 1 • - Photograph and schematic representation of shock tube test 
section for reflected shock waves. 


The information from the passage of the incident shock wave sensed by each 
transducer was recorded by a digital oscilloscope. The incident shock velocity, 
which was used to calculate the conditions behind the reflected shock wave, was 
obtained by using the time differential between the recorded signals of the trans- 
ducers located 19.23 cm and 3 mm from the end wall. These two transducers were 
separated by 18.93 cm. The measured velocities resulted in calculated, incident- 
shock-wave Mach numbers ranging from about 2 to 3. The time constant for the pres- 
sure measurements was 4 ps ± 0.5 jos for the entire series of tests. 

The shock tube was filled to a pressure of 40 mm of mercury with test gas for 
all the experiments. After approximately seven tests, the shock tube was disassem- 
bled and cleaned. This need for cleaning resulted partly from the debris of the 
diaphragms but mostly from the oxides of silicon powder formed from silane combus- 
tion. Once deposited on the shock tube walls, the powdered silicon oxides remained 
in place. 

Two gas mixtures with silane-oxygen equivalence ratios of 1.0 and 0.5 were used 
for the tests. The compositions of the test gas mixtures for silane-oxygen equiva- 
lence ratios of 1.0 and 0.5 are given in tables 1(a) and 1(b), respectively. 

These mixtures were obtained by diluting with nitrogen a commercially supplied 
mixture of 20 mole-percent silane and 80 mole-percent hydrogen. This SiH^-H 2 -N 2 
mixture was mixed with air in the shock tube to obtain the desired silane-oxygen 





stoichiometry in the test gas* This was implemented to prevent oxidation of the 
silane, which was observed to occur when premixing H^-SiH^-C^-^ at pressures above 
1 atm in the mixing vessel* When mixing was effected in the shock tube at pressures 
of 40 mm of mercury, no evidence of silane oxidation was noted (as would have been 
observed by the visible formation of a white silicon oxide deposit) • 


EXPERIMENTAL RESULTS AND DATA ANALYSIS 

An example of the pressure histories behind both the incident shock wave and the 
reflected shock wave is shown in figure 2. Figure 2(a) shows the pressure history 
measured by the transducer 19.23 cm from the end wall of the shock tube. As can be 
seen, the pressure quickly increases from the vacuum value to the pressure behind the 
incident shock wave (0.32 atm). Approximately 700 |js later, the shock wave, after 
reflection off the end wall, again passes the transducer and quickly raises the pres- 
sure level to approximately 1.25 atm. In figure 2(b), the pressure history 3 mm from 



(a) Transducer 19.23 cm from end wall. 



(b) Transducer 3 mm from end wall. 
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Figure 2.- Typical pressure histories for mixture II test gas after shock 

wave passage. Pressure = 1 .35 atm. 


the end wall may be observed. The first increase in pressure occurs when the shock 
wave initially encounters the transducer. For the end-on placed transducer, in all 
cases, the record was found to exhibit "ringing" or "overshoot." However , after 
passage of the shock wave, the pressure drops back toward the normal pressure level 
(*0.3 atm) behind an incident shock wave, as is labeled in the figure. After the 
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shock wave is reflected from the end wall of the shock tube (and travels the 3 mm), 
the pressure increases to the value of 1 ,25 atm behind a reflected shock wave. 
Approximately 200 \\s later (for the mixture illustrated), the mixture ignites. The 
ignition delay or induction time is measured from the point at which the pressure 
history exhibits its maximum slope after reflection of the shock wave, to the maximum 
slope reached after leveling off behind the reflected shock wave. 

The incident shock velocity needed in the calculation of the reflected shock 
wave conditions was obtained from the time difference between the transducers placed 
19.23 cm from the end wall and 3 mm from the end wall, as shown in figure 2. 

The measured ignition delay times % for the silane-hydrogen-air mixtures are 
shown in tables 11(a) and 11(b) along with the associated initial (unreacted) temper- 
atures behind the reflected shock wave. The mixtures containing silane required from 
24 to 752 jjs to ignite after passage of the reflected shock wave. The initial 
reflected shock temperatures were calculated from the incident shock speed with the 
assumption of no reaction behind the incident wave. These values ranged from approx- 
imately 800 K to 1050 K. The observed pressures just behind the reflected shocks 
were within 3 to 5 percent of those calculated, and the estimated uncertainty in the 
initial reflected shock temperature is no more than +10 K. The analytical model used 
to calculate the conditions behind the reflected shock wave is described in refer- 
ence 3. The variation of the experimentally measured ignition delay times over a 
range of temperatures is shown in figure 3 for a silane-oxygen equivalence ratio 



Figure 3.- Variation of ignition delay time with 
ture I. Si lane -oxygen equivalence ratio of 


reciprocal temperature for mix- 
1 .0; Pressure = 1 .25 atm. 
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of 1.0 and a pressure of 1.25 atm. The variation of the ignition delay times with 
temperature for a si lane -oxygen equivalence ratio of 0.5 and a pressure level 
of 1.35 atm is shown in figure 4. To further substantiate that silane enhances 
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Figure 4.- Variation of ignition delay time with reciprocal 

temperature for mixture II. Silane-oxygen equivalence ratio 
of 0.5; Pressure = 1 .35 atm. 


ignition, experiments with similar mixtures that did not contain silane were con- 
ducted over the same pressure and temperature range used with the si lane -hydrogen 
mixtures. Ignition did not occur during the available test times (»1000 \xs) • The 
lines in figures 3 and 4 represent least- squares curve fits to the data presented in 
tables 11(a) and 11(b). The curve fit shown in figure 3 resulted in the equation 

c -4 11200/T 

x = 5.71 x 10 e 


to describe the ignition delay time x (in microseconds) as a function of reflected 
shock temperature T (kelvin) for the mixture with a silane-oxygen equivalence ratio 
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of 1.0 and a test pressure of 1.25 atm. The least-squares curve fit of the data 
shown in figure 4, obtained with a silane-oxygen equivalence ratio of 0.5 and a test 
pressure of 1.35 atm, yielded the equation 


-4 12800/T 
T = 1 .18 x 10 e 


COMPARISON BETWEEN EXPERIMENTAL RESULTS AND KINETIC MODEL 

Beach et al. (ref. 2) proposed a mechanism for the combustion of silane which 
was assembled by analogy with the methane oxidation mechanism. The reactions and 
rate coefficients which made up the mechanism are listed in table III. Most of the 
rate coefficients for the silane reactions were estimated, and several were adjusted 
to fit the experimental results of Strater (ref. 7), who studied the combustion of 
silane in a flow reactor over the temperature range from 500 K to 700 K. 

Thermochemical properties were unavailable for some of the chemical species 
needed to develop a silane analog to a methane reaction mechanism. It was necessary, 
therefore, to estimate thermochemical property data for those species not in JANAF 
(ref. 4) or the open literature. As an example, the thermochemical data for the 
compounds SiH^ and SiH^ were estimated by analogy to the "families" of curves 
reported for hydrocarbons in reference 8. The assumption was made that silane is 
analogous to the methane family with species CH, CH 2 , CH^, and CH^, and the family of 
SiH, SiH^, SiH^, and SiH^ was generated. The data for SiH 2 and SiH^ were estimated 
to fit between the JANAF (ref. 4) data for SiH and SiH^, in a manner similar to that 
of reference 8. Other compounds, such as SiH^O and HSiO, were assumed to be present 
in the silane mechanism by analogy with methane. Thermochemical data for these com- 
pounds were estimated also. The analytical computer program described in reference 9 
was used to calculate the temperature and pressure history after passage of a 
reflected shock wave. To simulate the conditions after passage of a reflected shock 
wave, the computer program was operated in a constant volume (density) manner as 
described in reference 6. 

A comparison of the ignition delay measurements with the calculated ignition 
delay times for mixture I is shown in figure 5. Even though the mechanism predicts 
significant ignition enhancement over that for no silane (ref. 2), the measured 
enhancement is greater than predicted by the model. Considering the uncertainties in 
the silane kinetic mechanism and the fact that it was assembled and optimized by 
comparison with experimental results from a low temperature (500 K to 700 K) study of 
silane oxidation (ref. 7), the difference of a factor of 5 between the calculated and 
experimental results is reasonably good. However, the fact that a difference exists 
points out the need for further refinement of the reaction mechanism. In addition, 
better estimates of the thermochemical properties of the silicon-containing inter- 
mediates are needed to increase the reliability of the model. 
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l/T, K -1 

Figure 5,- Comparison of theoretically calculated ignition delay times with 
experimentally measured data for mixture I. 


CONCLUDING REMARKS 

* Ignition delay times for mixture I (2.00 percent silane, 8.00 percent hydrogen, 

4.01 percent oxygen, and 85.99 percent nitrogen) and mixture II (1.68 percent silane, 
6.72 percent hydrogen, 6.74 percent oxygen, and 84.86 percent nitrogen) have been 
< determined in a chemical kinetic shock tube with heating generated behind the 

reflected shock wave. The ignition delay times were measured at initial temperatures 
between approximately 800 K and 1050 K at pressures of 1.25 and 1.35 atm. Hie 
mixtures represented silane-oxygen equivalence ratios of 1 .0 and 0.5, respectively. 
Tests conducted with mixtures containing no silane failed to ignite over the same 
temperature-pressure regime. This definitely indicated enhancement of ignition 
through the use of silane. The mixtures containing silane required from 24 to 
752 ps to ignite after passage of the reflected shock wave. The measured ignition 
delay times were plotted against reciprocal temperature and curve -fitted by the least- 
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squares regression method. The equation for ignition delay time x (in micro- 
seconds) as a function of postshock temperature T (kelvin) for the silane-oxygen 
equivalence ratio of 1 .0 was 


x 


5.71 


-4 11200/T 
x 10 e 7 


For a silane-oxygen equivalence ratio of 0.5, the equation was 


x 


1 .18 x 


-4 

10 e 


12800/T 


Thermochemical properties and chemical kinetic coefficients were estimated and 
used in a computer model of the chemical system behind a reflected shock wave. The 
chemical kinetic package and estimated thermochemical properties, based on analogy to 
the methane system, yielded a comparison which was closer than expected to the exper- 
imentally observed values, especially considering that the estimated theoretical 
values were developed from a scarce amount of data. If the use of silane as a com- 
bustion initiator is to be investigated in greater detail, estimated thermochemical 
properties and reaction rate coefficients must be replaced with measured or calcu- 
lated values from spectrometric data. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 3, 1982 
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TABLE I.- COMPOSITION OF TEST GAS MIXTURES 


(a) Mixture I (silane-oxygen equivalence 
ratio of 1.0) 


Gas 

Mole-percent 

SiH 

2.00 

H 2 

8.00 

°o 

4.01 

z 

N 2 

85.99 


(b) Mixture II (silane-oxygen equivalence 
ratio of 0.5) 


Gas 

Mole-percent 

SiH 

1 .68 

H 2 

6.72 

°2 

6.74 

N 2 

84.86 
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TABLE II 


IGNITION DELAY TIMES 


(a) Mixture I (si lane -oxygen equivalence ratio of 1.0) 
[Pressure = 1.25 atm] 


Reflected shock 
temperature, T, K 

Ignition delay time, 
T, |0S 

1/T, K _1 

807 

635 

1.239 x 10 

815 

470 

1.220 

806 

635 

1 .241 

826 

433 

1 .211 

812 

460 

1 .232 

820 

443 

1.220 

845 

272 

1 .183 

882 

149 

1.134 

893 

137 

1 .120 

893 

170 

1.120 

893 

129 

1 .120 

876 

204 

1.142 

877 

142 

1 .127 

865 

204 

1 .1 56 

865 

187 

1 .156 

940 

94 

1 .064 

904 

128 

1 .106 

904 

128 

1 .106 

922 

24 

1 .105 

910 

135 

1.099 

866 

226 

1 .155 

887 

198 

1 .127 

860 

267 

1 .153 

870 

190 

1 .150 

835 

345 

1 .198 

850 

268 

1.177 

850 

300 

1 .177 

807 

689 

1 .239 

835 

400 

1 .190 

835 

452 

1.198 

846 

391 

1 .192 

81 1 

559 

1 .233 

807 

623 

1 .233 

81 1 

633 

1.233 

811 

573 

1 .232 

81 1 

628 

1.233 

801 

590 

1 .248 

807 

534 

1.239 

807 

558 

1 .239 

922 

93 

1 .085 

927 

100 

1 .079 

958 

74 

1.044 

962 

78 

1 .040 

976 

58 

1 .025 

983 

60 

1 .017 

990 

50 

1.010 

1046 

25 

.906 

1033 

25 

.968 

1046 

24 

.956 


v-3 



TABLE II.- Continued 


(b) Mixture II (silane-oxygen equivalence ratio of 0.5) 
[Pressure = 1.35 atm] 



1/T, K" 1 


1.075 x 10 
.969 
1 .022 
1 .202 
1 .193 
1.176 
1 .030 
1 .078 
1 .035 
.982 
1 .042 
.976 
1 .099 
1 .091 
1 .055 
1 .050 
1 .163 
1 .099 
1 .142 
1.149 
1 .176 
1 .201 
1 .163 
1 .190 
1 .222 
1 .217 
1 .149 
1.105 
1 .163 
1.163 
1 .105 
1 .078 


TABLE III 


SILANE-HYDROGEN REACTION MECHANISM 


[From ref. 2] 


Reaction 

*A 

*n 


+ M + SiH 4 -*■ SiH 3 + H + M 

1 .00 X 10 18 

0 

68 051 

SiH 4 + H -»■ SiH 3 + H 2 

2.50 x 10 15 

0 

3 700 

SiH 4 + 0 -> SiH 3 + OH 

8.50 x 10 13 

0 

3 850 

SiH 4 + OH SiH 3 + H 2 0 

5.70 x 10 13 

0 

5 000 

SiH 3 + 0 2 ■> SiH 2 0 + OH 

5.70 x 10 13 

0 

7 000 

SiH 3 + 0 -► HSiO t H 2 

1.00 x 10 15 

0 

0 

M + SiH 2 0 -> HSiO t H + M 

5.00 x 10 12 

.736 

30 000 

SiH 2 0 + H •> HSiO + H 2 

2.00 x io 15 

0 

800 

SiH 2 0 + 0 ■> HSiO + OH 

to 

. 

o 

o 

X 

o 

in 

0 

3 500 

SiH 2 0 + OH -*• HSiO + H 2 0 

2.00 x io 15 

0 

0 

M + HSiO -> SiO + H + M 

1.00 x 10 15 

0 

23 000 

HSiO + H ->■ SiO + H 2 

1 .00 x to 15 

0 

2 300 

HSiO + 0 -> SiO + OH 

1 .00 x 10 15 

0 

2 400 

HSiO + OH -*• SiO + H 2 0 

1 .00 x 10 15 

0 

1 400 

Si0 2 + OH ->■ HSiO + 0 2 

7.20 x io 10 

.5 

80 000 

SiO + H0 2 -> HSiO + 0 2 

. 

'O 

o 

X 

o 

.5 

4 900 

SiO + OH -> Si0 2 + H 

4.00 X 10 15 

0 

4 300 

SiO + 0 2 ■> Si0 2 + 0 

1 .00 X io 14 

0 

5 200 


*The parameters A, n, and E^ refer to the rate coeffi- 
cient equation k = AT n exp(-E^/RT). The rate coefficient units 
are s’* 1 for unimolecular reactions, cm 3 /mole-s for bimolecular 
reactions, and cm 6 /mole 2 -s for termolecular reactions. The 
symbol E^ is the activation energy (cal/mole) , R is the gas 
constant (1 .987 cal/mole-K) , and n is a temperature exponent 
factor obtained from empirical correlation. 

+ M is a third-body molecule. 



TABLE III,- Concluded 


[From ref* 2] 


Reaction 

*A 

*n 


SiO + H0 2 Si0 2 + OH 

1.00 X 10 14 

0 

2 000 

SiO + 0 + + M ->■ Si0 2 + M 

6.00 x 10 16 

0 

0 

H + 0 2 + M > H0 2 + M 

2.10 x 10 18 

-1 

0 

H 2 + OH •> H 2 0 + H 

2.19 x 10 13 

0 

5 150 

H + 0 2 ■* OH + 0 

2.20 x 10 14 

0 

16 800 

0 + H 2 -*• OH + H 

1 .80 x 10 10 

1 

8 900 

H 2 + H0 2 -*• H 2 0 2 + H 

3.00 x 10 11 

0 

18 680 

M + H 2 0 2 -* OH + OH + M 

1.20 x io 17 

0 

43 500 

H + H0 2 OH + OH 

2.00 x io 14 

0 

1 800 

H + OH + M •* H 2 0 + M 

8.40 x io 21 

-2 

0 

o + o + m>o 2 + m 

1.38 x io 18 

-1 

340 

0 + H 2 0 •> OH + OH 

6.80 x io 13 

0 

18 360 

H 2 + ° 2 °H + OH 

1 .00 x 10 13 

0 

43 000 

OH + H0 2 ■* H 2 0 + 0 2 

5.00 x io 13 

0 

1 000 

O + H0 2 + OH + 0 2 

5.00 x io 13 

0 

1 000 

H°2 + H °2 ■” H 2°2 + °2 

2.00 x io 12 

0 

0 

H + H0 2 H 2 + 0 2 

2.50 x 10 13 

0 

700 

h+h + m->-h 2 + m 

6.40 x 10 17 

-1 

0 


*The parameters A, n, and ^ refer to the rate coeffi- 
cient equation k = AT n exp(-E R /RT). The rate coefficient units 
are s” 1 for unimolecular reactions, cm 3 /mole-s for bimolecular 
reactions, and cm 6 /mole 2 -s for termolecular reactions. The 
symbol E^ is the activation energy (cal/mole), R is the gas 
constant (1.987 cal/mole-K), and n is a temperature exponent 
factor obtained from empirical correlation, 
is a third-body molecule. 



1. Report No. 2. Government Accession No. 

NASA TP-2114 

3. Recipient's Catalog No. 

4. Title and Subtitle 

IGNITION OF SiH 4 -H 2 -0 2 -N 2 BEHIND REFLECTED 
SHOCK WAVES 

5. Report Date 

February 1983 

6. Performing Organization Code 
505-43-83-03 

7. Author(s) 

Allen G. McLain, Casimir J. Jachimowski, and 
R. Clayton Rogers 

8. Performing Organization Report No. 
L- 1 5534 

10. Work Unit No. 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Paper 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 


The ignition of mixtures of silane, hydrogen, and oxygen diluted with nitrogen was 
studied in a chemical kinetic shock tube. Ignition delay time measurements were made 
behind reflected shock waves at pressures of 1,25 atm and 1.35 atm and temperatures 
ranging from approximately 800 K to 1050 K for two mixtures with silane-oxygen equiv- 
alence ratios of 1.0 and 0.5, respectively. Enhancement of the ignition by use of 
silane was observed experimentally by the failure of a similar mixture containing 
only hydrogen, oxygen, and nitrogen to ignite during the maximum test time available 
for these test conditions. 


17. Key Words (Suggested by Author(s)) 


18. Distribution Statement 


Ignition delay times 
Silane ignition 
Shock ignition 


Unclassified - Unlimited 


Subject Category 07 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

Unclassified 

Unclassified 

15 


22. Price 
A0 2 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1983 



National Aeronautics and 
Space Administration 


THIRD-CLASS BULK RATE 


Washington, D.C. 
20546 

Official Business 

Penalty for Private Use, $300 


1 iU ’ A * 330214 

OtFl Or 1.3- A In' FOwCE 
Ah »Vh AEONS LA.30iMro.2y 
ATO* I'cCH.VI CAL LIBaAa/ 
<Ii 21L.AU ARi N'.'i 37117 


330 y 3 303 


(S-JL) 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 








POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 



